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ABSTRACT 



i 

A toy model for gamma-ray burst-supernovae (GRB-SN) is discussed by considering the 
coexistence of baryon poor outflows from black holes (BHs) and a powerful spin-connection 
to surrounding disk, giving rise to consistent calorimetry as described by van Putten (2001a) 
in a variant of the Blandford-Znajek process (BZ, 1977). In this model the half-opening angle 
of the magnetic flux tube on the horizon is determined by the mapping relation between the 
ON I angular coordinate on the BH horizon and the radial coordinate on the surrounding accretion 
disk. The GRB is powered by the baryon poor outflows in the BZ process, and the associated 
SN is powered by very small fraction of the spin energy transferred from the BH to the disk 
in the magnetic coupling (MC) process. The timescale of the GRB is fitted by the duration 
of the open magnetic flux on the horizon. It turns out that the data of several GRB-SNe 
are well fitted with our model. 



O 

S3 

Subject headings: black hole physics-accretion disks-gamma-ray bursts-supernovae 

CO 

1. INTRODUCTION 

X 

As is well known, the Blandford-Znajek (BZ) process is an effective mechanism for powering jets 
from quasars and active galactic nuclei (Blandford & Znajek 1977; Rees 1984). In the BZ process energy 
and angular momentum are extracted from a rotating balck hole (BH), and transferred to the remote 
astrophysics load by the open magnetic field lines. Recently, much attention has been paid on the 
issue of the long duration gamma-ray bursts (GRBs) powered by the BZ process (Paczynski 1993, 1998; 
Meszaros & Rees 1997). A detailed model of GRBs invoking the BZ process is given by Lee, Wijers, 
Sz Brown 2000, hereafter LOO), and they concluded that a fast-spinning BH with strong magnetic field 
~ 10 15 G can provide energy ~ 10 53 ergs to power a GRB within 1000 seconds. 

Recently the observations and theoretical considerations have linked long duration GRB with ultra- 
bright type Ib/c supernovae (Galama et al. 1998; Bloom et al. 1999; Galama et al. 2000). The first 
candidate was provided by SN1998bw and GRB980425, and the recent HETE-II burst GRB030329 has 
greatly enhanced the confidence in this association (Stanek et al. 2003; Hjorth et al. 2003). 
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Not long ago Brown et al. (2000, hereafter BOO) worked out a specific scenario for GRB-SN 
connection. They argued that the GRB is powered by the BZ process, and the SN is powered by the 
energy dissipated into the disk through closed magnetic field lines coupling the disk with the BH. The 
latter energy mechanism is referred to as the magnetic coupling (MC) process, which is regarded as one 
of the variants of the BZ process (Blandford 1999; Putten 1999; Li 2000, 2002; Wang, Xiao & Lei 2002, 
hereafter W02). It is shown in BOO that about 10 53 ergs are available to power both GRBs and SNe. 
However, they failed to distinguish the fractions of the energy for these two objects. 

More recently, van Putten et al. (van Putten 2001a; van Putten & Levinson 2003; van Putten 
et al. 2004) worked out a poloidal topology for the open and closed magnetic field lines, where the 
separatrix on the horizon is defined by a finite half-opening angle. The duration of GRB is set by the 
lifetime of rapid spin of the BH. It is found that GRB and SN are powered by a small fraction of the 
BH spin energy. This result is consistent with the observations, i.e., the duration of GRB with tens of 
seconds, the true GRB energies distributed around 5 x 10 50 ergs (Frail et al. 2002), and the aspherical 
SNe kinetic energies with 2 x 10 51 ergs ( Hoflich et al. 1999). 

In this paper we propose a toy model for GRB-SN by considering coexistence of the BZ and MC 
processes (CEBZMC), where the configuration of the magnetic field is based the works of Wang et al. 
and van Putten et al. (W02; Wang et al. 2003, hereafter W03; van Putten 2001a; van Putten et al. 
2004). The transferred energy and the duration for GRBs are calculated in the evolving process of 
the half-opening angle of the magnetic flux on the BH horizon, while the transferred energy for SNe 
is calculated in the MC process. Thus the association of GRBs and SNe is explained reasonably and 
consistently. 

This paper is organized as follows. In § 2 the configuration of the magnetic field is described with 
the half-opening angle 6bz of the magnetic flux tube on the horizon, and we find that 9bz can evolve to 
zero in some range of the power-law index n of the magnetic field on the disk. In § 3 the toy model for 
GRB-SN based on BH evolution is described in the parameter space consisting of the BH spin a* and the 
power-law index n. It is shown that the half-opening angle 9bz plays an important role of determining 
(i) the powers for GRBs and SNe, (ii) the duration of GRB. In § 4 we calculate the duration of GRB 
and the energies of GRB-SN, and compare with the other models of GRBs invoking the BZ process. 
It is found the duration of GRB obtained in our model is shorter than that of the other models. It 
turns out that this model is in excellent agreement with the observations, owing to considering the MC 
effects, while too much energy than the needed energy of GRB is produced in the other models without 
the MC effects. Finally, in § 5, we summarize our main results and discuss some problems concerning 
our model. Throughout this paper, the geometric units G = c = 1 are used. 



2. CONFIGURATION OF THE MAGNETIC FIELD 

Recently, we proposed a model of CEBZMC, in which the remote astrophysical load in the BZ 
process and the disk load in the MC process are connected with a rotating BH by open and closed 
magnetic field lines, respectively (W02; W03). The poloidal configuration of the magnetic field is shown 
in Figure 1, which is adapted from van Putten (2001a). 

In Figure 1 the angle Qbz is the half-opening angle of the open magnetic flux tube, indicating the 
angular boundary between open and closed field lines on the horizon. This angle has been discussed 
by van Putten and his collaborators, who related the half-opening angle to the curvature in poloidal 
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topology of the inner torus magnetosphere (van Putten & Levinson 2003; van Putten et al. 2004). In 
W03 the angle 6 B z is determined based on the following assumptions. 

(1) The theory of a stationary, axisymmetric magnetosphere formulated in the work of MacDonald 
& Thorne (1982) is applicable not only to the BZ process but also to the MC process. The magnetosphere 
is assumed to be force-free outside the BH and disk. 

(2) The disk is both stable and perfectly conducting, and the closed magnetic field lines are frozen 
in the disk. The disk is thin and Keplerian, and it lies in the equatorial plane of the BH, with inner 
boundary being at the marginally stable orbit. 

(3) The poloidal magnetic field is assumed to be constant on the horizon and to vary as a power 
law on the disk as follows (Blandford 1976), 

b d oc r n , (i) 

where Bd is the magnetic field on the disk, the parameter n is the power-law index, and £ = r/r ms is 
the radial coordinate on the disk, which is defined in terms of the radius r ms = Mx 2 ms °f the marginally 
stable orbit (Novikov & Thorne 1973). 

(4) The magnetic flux connecting the BH with its surrounding disk takes precedence over that 
connecting the BH with the remote load. 

Assumption (4) is proposed based on two reasons: (i) the magnetic field on the horizon is brought 
and held by the surrounding magnetized disk, and (ii) the disk is much nearer to the BH than the 
remote load. 

A mapping relation between the parameter £ and the angular 9 is derived in W03 based on the 
conservation of magnetic flux of the closed field lines. 

cos9 = J* G (a*; £,n)dS, (2) 

where a* = J/M 2 is the BH spin defined by the BH mass M and angular momentum J, and the function 
G(a*] £, n) is given by 



G(a*;£,n) = . (3) 

V(! + a lX~ m \ + 2alxJ s ) (1 - 2x- 2 s £" 1 + aSx^" 2 ) 

Based on the assumption (4) the angle Obz can be determined by taking £ = oo for the highest 
closed field lines as follows, 



cos# 



BZ 



/oo 
G(o,;£n)de. (4) 



The curves of Obz versus a* with different values of power-law index n are shown in Figure 2. 
Inspecting Figure 2, we have the following results. 
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(1) The half-opening angle 9bz always increases monotonically with the increasing n. 

(2) For some values of n, such as n— 3.5, 4, 4.5, 5 and 5.493, the angle 9bz can evolve to zero with 
the decreasing a*. 

The second point implies that the open magnetic flux tube will be shut off when the BH spin 
decreases to the critical value a^ RB corresponding to 9bz = 0. The lifetime of the half-opening angle 
9bz is defined as the evolution time of BH from the initial spin a*(0) to a^ RB . Obviously it is less than 
the lifetime of BH spin, which is the evolution time of BH from a*(0) to zero. By using equation (4) we 
find that a^ RB exists only for a specific value range of the power-law index n, i.e., 3.003 < n < 5.493. 
The half-opening angle Qbz remains zero for n < 3.003 in the whole evolving process, and it will never 
evolve to zero for n > 5.493. 



3. A TOY MODEL FOR GRB-SN BASED ON BH EVOLUTION IN CEBZMC 

Considering that the angular momentum is transferred from a rapidly rotating BH to the disk, on 
which a positive torque exerts, we think the accretion onto the BH might be probably halted. This 
state is essentially the suspended accretion state proposed by van Putten & Ostriker (2001), and the 
evolution of the BH is governed by the BZ and the MC processes. 

The expressions for the BZ and MC powers and torques are derived in W02 by using an improved 
equivalent circuit for the BH magnetosphere based on the work of MacDonald and Thorne (1982). 
Considering the angular boundary 9bz, we express the BZ and MC powers and torques as follows. 

P ~P IP ~2, 2 f° BZ k (I -k) sin 3 Odd 

Pbz = Pbz/Po - 2a, / 2 _ (1 _ g) ^ (5) 



o 



[1 - k) sin 3 9d9 
2- (l-q) sin 2 9' 



f B z = T BZ / To = 4a, (! + ?)/ ^ * } ( 6 ) 



2 r /2 (3(1 -f3) sin 3 9d9 
le BZ 2 -{l-q) sin 2 # 



Pmc = Pmc/Po = 2a 2 / ^_ ^ , (7) 



t — rp I rp A h , x r' 2 (1-13) sin 3 9d9 

Tmc = T MC /T = 4a, 1 + q) / - — r-^, 

1 Je BZ 2 - (1 - q) sin 9 

where we have q = ^/l — a 2 and 

P = (B H ) M 2 « 6.59 x 10 50 x {^f erg ■ 

To = (B 2 H ) M 3 « 3.26 x 10 45 (^)% • cm 2 • S " 



(8) 



(9) 



In equation (9) B>h is the magnetic field on the BH horizon. The parameters k and (3 are the ratios 
of the angular velocities of the open and closed magnetic field lines to that of the BH, respectively. 
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Usually, k = 0.5 is taken for the optimal BZ power. Since the closed field lines are assumed to be frozen 
in the disk, the ratio (3 is related to a* and £ by 

p = n D /n H = r 2(1 + ( ? — T , (io) 



(v^Xm S ) + a* 

where and Op are the angular velocities of the BH and the disk, respectively. 

Since 9 BZ in equations (5) and (7) is determined by a* and n, we have the curves of the powers 
Pbz and Pmc versus a* for the given values of n as shown in Figure 3. 

Inspecting Figure 3, we have the following results. 

(1) The MC power Pmc is always greater than the BZ power P BZ for a wide value range of the 
power-law index n, provided that the BH spin a* is not very small. 

(2) For the given values of a* the BZ power P BZ increases with the increasing n, while the MC 
power Pmc decreases with it. 

(3) The power P M c = holds at some small values of a* for the given values of the index n. 

From equation (5) we find Pbz = when 6bz = 0. Therefore the effective time of the BZ process 
can be set by the lifetime of Qbz, i.e., the evolution time of BH from a*(0) to a^ RB for 3.003 < n < 5.493. 
As argued in our previous work (Wang, Lei & Ma 2003), we have Pmc — when the energy transported 
from the BH to the disk is equal to that from the disk to the BH. Henceforth the corresponding critical 
BH spin is denoted by af N . We have Pmc > for a* > af N , and Pmc < for a* < af N , indicating 
that the energy transferred from the BH to the disk dominates that transferred in the inverse direction. 

Based on the conservation laws of energy and angular momentum we have the following evolution 
equations of the rotating BH, 

dM/dt=-(P BZ + PMc), (11) 



dJ/dt = -(Tbz + Tmc)- (12) 
Incorporating equations (11) and (12), we have the evolution equation for the BH spin expressed 

by 

dajdt = -M~ 2 (T BZ + T MC ) + 2M- 1 a*(P BZ + Pmc) = B 2 H MA(a*, n) (13) 
where the function A(a*,n) is 

where 



P-mag ~ PbZ + PmCi T mag — T BZ + T MC - 



(15) 



The sign of A (a*, n) determines whether a* decreases or increases, and we have the curves of A(a*,n) 
versus a* for different values of n as shown in Figure 4. 

From Figure 4, we find the following evolution characteristics of the BH spin. 

(1) The function A(a*,n) > for a* < a eq , and A(a*,n) < for a* > al q , and the BH spin always 
evolve to the equilibrium spin a eq . 

(2) Both the equilibrium spin a eq and the decreasing rate of the BH spin are sensitive to the power- 
law index n. The greater the index n is, the greater is the equilibrium spin a% q , and the more slowly the 
BH spin decreases to a eq . 

(3) The features of our model for GRB-SN can be described by three functions: A(a*, n), Pbz(<i*, n) 
and P MC (a*,n). 

Setting A(a*,n) = 0, PBz(a*,n) = and PMc( a *, n ) = 0, we have the three characteristic curves 
in the parameter space as shown in Figure 5. 

The parameter space is divided into seven regions, IA, IB, II, IIIA, IIIB, IVA and IVB, where 
the points M and N are the intersections of the curve PBz(a*,n) = with the curves PMc( a *, n ) — 
and A(a*,n) = 0, respectively. IA and IB is divided by the dashed line starting from the point M. By 
using the equations (5), (7) and (14) we obtain the coordinates of points M and N, i.e., M(0.253, 
5.067) and 7V(0.222, 5.122), respectively. In Figure 5 each black dot with an arrowhead in these 
sub-regions is referred to as a representative point (RP), which represents one evolution state of the 
BH as shown in Table 1. 

It is found from Figure 5 and Table 1 that GRB-SN only occurs in region I and I corresponding 
to CEBZMC. Based on the observations that GRBs need less energy than SNe (BOO), we think that 
the evolving path of the RP in region IA is reasonable, which corresponds to 3.003 < n < 5.067, and 
the RP attains PBz(a*,n) = with PMc( a *, n ) > in braking the spinning BH. 

Based on the above discussion we describe the toy model for GRB-SN in the parameter space as 
follows. 

(1) The fast-spinning BH can power GRBs and SNe via the BZ and MC processes, respectively, 
which is represented by the RP in region I A. The GRB will be shut off, when the RP reaches eventually 
the curve Pez(o*,n) = 0. 

(2) Energy is transferred from the spinning BH to the disk via the MC process in the evolving path 
from region IA to region II. The SN will be shut off with no further energy transferred into the disk, 
if the RP arrives at the curve P MC (a*,n) = 0. 

(3) Since enormous energy is deposited in the disk due to the MC process, the disk might be de- 
stroyed eventually, leaving the BH alone without the disk. This outcome probably arises from explosion 
of SNe (BOO). 
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4. ENERGIES AND TIME SCALES FOR GRB-SN 

In the evolving path of the RP going through regions IA, II and IIIA we have af RB > af N > a% q . 
Based on the above discussion on correlation of the BH evolution with the association of GRB-SN , the 
energy Ebz and Esn extracted in the BZ and MC processes are given respectively by 



r Pbz , 17Q in54 (M(o)\ n* MP BZ 

E B z = / n — nrrda* = 1.79 x W^ergs x ; / — = = da*, (16) 



E MC = f " Tptrrrda. = 1.79 x lO^ergs x ( f " , MPmC . da*. (17) 
Jajo) [da*/dt) \ M & J J , 0) -T +2a*P maq 



The true energy for GRBs, E 7 , and the energy for SNe, E S n, are related respectively to E B z and 
E M c by 

i? 7 = e^Ebz, (18) 

-Esiv = £snEmc, (19) 

where e 7 and egjv denote the efficiencies of converting E B z and E MC into _E 7 and E SN , respectively. 
Following van Putten et al. (2004), we take e 1 = 0.15 in calculations. The duration of GRB, t GRB , is 
defined as the lifetime of the half-opening angle 9bz, which is exactly equal to the time for the BH spin 
evolving from a*(0) to af RB , i.e., 



tcRB = 2.7 x 10 3 s x — — — g_ / = rfa„ 

I ^ y \M(0))J a4o) -T mag + 2a*P mag 

where M is the ratio of the BH mass to its initial value M(0), and it is calculated by 



(20) 



M f a * da* . , 

M = — — = exp / j—^ . (21) 



M(0) 



a * (°) T mag / P ma g — 2a, 



Equation (21) can be derived by using equations (11) and (13). Obviously, t GR B is shorter than the 
lifetime of the BH spin evolving from a*(0) to a* = 0. 

Considering that spin-energy of the BHs produced in core-collapse is around 50% of maximum or 
less in centered nucleation (van Putten 2004), we take the initial BH spin as a*(0) = 0.9 in equations 
(16), (17) and (21). The value range 3.534 < n < 5.067 is taken in the evolving path of the RP in 
region IA corresponding to GRB-SN. In addition, M(0) = 7M and B H = 10 15 G are assumed, and the 
cutoff of T GRB is taken as T 90 in calculations, which is the duration for 90% of the total BZ energy to 
be extracted (Lee & Kim 2002). 
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As is well know, the rotational energy E rot of a Kerr BH is expressed by (Thorne, Price & Macdonald 
1986) 

E rot = /(a*)M(0), (22) 

where 

/(«*) = i - v^ + aA 3 ^!)- (23) 

Substituting a*(0) = 0.9 and M(0) = 7M into equation (22), we have the total rotational energy of 
the BH, 

E rot « 1.9 x I0 54 ergs. (24) 

To compare rotational energy of the BH extracted in the BZ and MC processes, we calculate the 
ratios of Ebz and Emc to E rot , i.e., 

7]bz = E B z/E rot , 7]mc = E M c/E rot . (25) 

The curves of r/sz and t/mc versus n are shown in Figure 6. 

From Figure 6 we find that -Emc* is significantly greater than Ebz, and the greater the indexn, 
the more energy extracted by BZ process. These results can be also obtained by inspecting Figures 3 
and 5. The ratio f]MC varies from 25.8% to 33.8% with the ratio t]bz varying from zero to 3.7% for 
3.534 < n < 5.067. Combining equations (24) and (25) with the above ratios, we obtain Emc varying 
from 4.929 x 10 53 ergs to 6.459 x 10 53 ergs with E B z = 7.003 x W 52 ergs. And the energy for GRB is 
inferred by using e 7 = 0.15. 

Calculations of SN1998bw with aspherical geometry show that the needed kinetic energy is about 
2 x 10 51 ergs (Hoflich et al. 1999), and this results in the ratio e S n varying from 0.003 to 0.004 for E MC 
varying from 6.459 x 10 53 ergs to 4.929 x 10 53 ergs in equation (19). In this paper we take the average 
value, Esn ~ 0.0035. 

It has been shown that most fraction of the rotational energy of a Kerr BH is emitted in unseen 
channels, such as in gravitational radiation and MeV-neutrino emissions (van Putten 2001a, 2001b; 
van Putten and Levinson 2002). The above efficiency factors obtained in our model, e 1 = 0.15 and 
e S N ~ 0.0035, imply that only very small fraction of E rot is converted into the energy for GRBs and 
SNe. These results are consistent with those obtained by van Putten and his collaborators. 

By using equation (20) the duration of GRBs, T 90 , versus the parameter n is shown in Figure 7. 

We find that the estimated duration of GRB is tens of seconds, which is consistent with the 
observations (Kouveliotous et al. 1993). For n = 4 and B H = 1 x 10 15 G in the evolving path of the RP 
from a*(0) = 0.9 to af N = 0.235, we obtain the duration of GRB T 90 = 21s, and the energies for GRB 
and SN are E 1 = 5 x 10 50 ergs(-^) and Esn = 2 x 10 51 ergs( £ Q^ 5 ), respectively. These results are in 
excellent agreement with the observed durations of tens of seconds (Kouveliotous et al. 1993), energies 
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Ej = 5 x 10 50 ergs in gamma rays (Frail et al. 2002), and inferred kinetic energy Esn = 2 x 10 51 ergs 
in SN1998bw with aspherical geometry (HofTich et al. 1999). 

For several GRB-SNe, the observed energy E 1 and duration T 90 can be fitted by adjusting the 
parameters n and Bh, and the energy Esn can be predicted as shown in Table 2. 

In the previous model for GRBs powered by the BZ process, the duration of GRB is estimated 
either in the case that the BH is spun down to zero or in the case that the whole disk is plunged into 
the BH (LOO; Lee & Kim 2000, 2002; Wang, Lei & Xiao 2002). However, the MC effects have not been 
taken into account in these models. In this model the duration of GRB is estimated by the lifetime of 
the half-opening angle on the BH horizon based on CEBZMC. It turns out that both the BZ and MC 
processes play very important roles for GRN-SN in this model. 

The quantities, E 7 and T 90 , calculated in the three different models for GRBs are listed in Table 
3, where the abbreviations CEBZMC, BZO and BZACC represent this model, the model invoking 
the BZ process only, the model invoking the BZ process with disk accretion, respectively. 

It is found in Table 3 that the energy and duration time obtained in this model are in excellent 
agreement with the observed durations of tens of seconds and energies E 1 m 5 x 10 50 ergs. However the 
energies obtained in the models BZO and BZACC seem too much to fit the observations. 



5. DISCUSSION 

In this paper we discuss a toy model for GRB-SN, where the energy for GRBs and SNe are powered 
by the BZ and MC processes, respectively, and the duration for GRB is estimated by the lifetime of 
the half-opening angle 6 B z- The energy extracted in the BZ process for GRBs is much less than that 
extracted in the MC process for SNe. This result is consistent with the observations: the true energy 
of GRB E 1 = 5 x 10 50 ergs and the aspherical SNe of kinetic energy Esn = 2 x 10 51 ergs. For a set of 
GRBs the observed true energy E 1 and duration T 90 can be well fitted by taking adequate values of the 
power-law index n and magnetic field Bh- However, there are still several issues related to this model. 

(i) In this paper we take the evolving path of the RP in the region IA to describe the association of 
GRB-SN as shown in Figure 5. In this case the GRB will terminate before the SN event stops. In fact 
the evolving path in region IB can be also used to describe the association of GRB-SN, and it implies 
that the RP will arrive at the characteristic curve of Pmc — ahead of reaching the characteristic curve 
of Pbz = 0. However, we think that the GRB will still terminate before the SN event stops, because 
the disk will be destroyed during the explosion of the SN, and the BZ process can not work without the 
magnetic field supported by its surrounding disk. 

(ii) In order to highlight the effects of the BZ and MC processes in powering GRB-SN, we neglect 
the effects of the disk accretion in this model. In fact, disk accretion plays a very important role 
in BH evolution. The slowing down of the BH spin will be delayed significantly since the accreting 
matter delivers remarkable amount of energy and angular momentum to the BH. Strictly speaking, disk 
accretion cannot be halted steadily by the MC process. The MC effects on disk accretion is a very 
complicated, and it should be treated dynamically. 

(iii) According to Kruskal-Shafranov criteria (Kadomtsev 1966), the screw instability will occur if 
the toroidal magnetic field becomes so strong that the magnetic field lines turns around itself about once. 
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In our latest work (Wang et al. 2004), we argued that the state of CEBZMC always accompanies the 
screw instability. In this model the half-opening angle 6bz is related to the infinite radial parameter £ by 
equation (4). We expect to derive a greater 6bz by considering the restriction of the screw instability to 
the MC region. Thus the ratios i] GRB and i] SN will be changed. We shall improve this model in future. 
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Table 1. Characteristic of BH evolution in each region of the parameter space 



Region 


da t 1 dt 


Pbz 


Pmc 


RP displacement 


Possible events 


IA 


< 


> 


> 


Towards the left 


GRB-SN 


IB 


< 


> 


> 


Towards the left 


GRB-SN 


II 


< 




> 


Towards the left 


SN 


IIIA 


< 




< 


Towards the left 




IIIB 


< 


> 


< 


Towards the left 


GRB 


IVA 


> 




< 


Towards the right 




IVB 


> 


> 


< 


Towards the right 


GRB 



Table 2. Four GRBs of the given true energy E 1 and T 90 s fitted with the different power-law index n 

and the predicted E SN 



GRB a 


E-y(W 51 ergs) h 


r 90 (s) 


n 


B„(10 15 G) 


E S N C 


970508 
990712 
991208 
991216 


0.234 
0.445 
0.455 
0.695 


15 d 
30 c 
39.84 f 
7.51 f 


3.885 
3.975 
3.985 
4.058 


0.97 
0.81 
0.75 
1.80 


1.947xl0 51 er 9 s( g ^) 
1.989xl0 51 er S s(^g) 
1.995xl0 51 er 9 s(^) 
2.032xl0 51 er 9 s(^) 



a The sample of GRB-SN takes from the sample listed in Table 1 of Dar 
(2004) 

b The true energy E 7 of GRB refers to the sample listed in Table 1 of Frail 
et al (2002) 

c The predicted energy of SNe based on our model 

d The duration of GRB970508 refers to Costa et al (1997) 

e The duration of GRB990712 refers to Heise et al. (1999) 

f The durations of GRB991208 and GRB991216 refer to the sample listed 
in Table 2 of Lee & Kim (2002) 



Table 3. Energy and duration obtained in three models for GRBs with a*(0) = 0.9, M H (0) = 7M , 

B H = 10 15 G and the initial disk mass M D (0) = 3M 



Model 


Ej (ergs) 


790 (s) 


BZO 


1.225 x 10 53 


348 


BZACC 


1.734 x 10 53 


397 


CEBZMC 


< 1.050 x 10 52 


< 115 



-13- 



GRB 




Black Hole Accretion Disk 



Fig. 1. — Poloidal magnetic field connecting a rotating BH with remote astrophysical load and a 
surrounding disk 
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Fig. 3. — Curves of Pgz (dotted lines) and P M c(solid lines) vs. a* with n=3.5, 4.5 and 5.493 for 
< a* < 1. 
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1 




Fig. 4. — Curves of A(a*,n) vs. a* for ra=3.5, 4, 4.5 and 5.493. 
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Fig. 5. — Parameter space of BH evolution in CEBZMC. 
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Fig. 6. — Curves of ^^(dashed line) and r] MC (so\id line) vs. n for 3.534 < n < 5.067. 
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